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ABSTRACT 

We present the result of 6 years monitoring of SgrA*, radio source associated to the 
supermassive black hole at the centre of the Milky Way. Single dish observations were 
performed with the Itapetinga radio telescope at 7 mm, and the contribution of the 
SgrA complex that surrounds SgrA* was subtracted and used as instantaneous cali- 
brator. The observations were alternated every 10 min with those of the HII region 
SrgB2, which was also used as a calibrator. The reliability of the detections was tested 
comparing them with simultaneous observations using interferometric techniques. Dur- 
ing the observing period we detected a continuous increase in the SgrA* flux density 
starting in 2008, as well as variability in timescales of days and strong intraday fluc- 
tuations. We investigated if the continuous increase in flux density is compatible with 
free-free emission from the tail of the disrupted compact cloud that is falling towards 
SgrA* and concluded that the increase is most probably intrinsic to SgrA*. Statistical 
analysis of the light curve using Stellingwerf and Structure Function methods revealed 
the existence of two minima, 156 ± 10 and 220 ± 10 days. The same statistical tests 
applied to a simulated light curve constructed from two quadratic sinusoidal functions 
superimposed to random variability reproduced very well the results obtained with the 
real light curve, if the periods were 57 and 156 days. Moreover, when a daily sampling 
was used in the simulated light curve, it was possible to reproduce the 2.3 GHz struc- 
ture function obtained by Falcke in 1999, which revealed the 57 days period, while the 
106 periodicity found by Zhao et al in 2001 could be a resonance of this period. 
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1 INTRODUCTION 



SgrA* is a compact radio source, with radius smaller than 1 
AU, which is associated to the 4x lO^M© supermassive black 
hole localized in the centre o f the Milky Way ([Eckart Sz Gen- 



zel|1996';'Menten et al.[T997l|Schodel et al.|2002||Eisenhauer 

et al. 2003; Gilles sen et al.|2009 ). Its mass and distance make 
the event horizon of SgrA* the largest angular size of all 
known black hole candidates. It is a sub-Eddington source, 
with bolometric luminosity 10~^° times its Eddington lumi- 
nosity (see Melia & Falcke (2001), for a review). Recently, 
Gillessen et al.| ( |2012t detected a gas cloud that is approach- 



ing SgrA* in a highly eccentric Keplerian orbit. The mass of 
the cloud has been estimated as three earth masses and will 
reach the black hole in the middle of 2013, during pericen- 
tre passage. The cloud is being elongated along the orbital 
path as a consequence of tidal disruption and will finally 
disintegrate into small clumps. 

SgrA* has been observed since the discovery of the radio 



source by Balick & Brown (1974) and its emission presents 



variability at all wavelengths, with time scales that can be 



just a few hours in the infrared (Genzel et al. 2003 Ghez 



et al. 



2004 



Meyer et al. 



( Baganoff et al. 20011 2003 ) , or days and months at radio fre- 



2008 



Do et al. 2009 ) and X-rays 



quencies ( Zhao et al.|1992' Duschl Lesch|1994 Herrnstein 



et al.||2004 ). Intra day variability (IDV) have been detected 
at radio wavelengths: 0.35, 0.7 and 1.3 cm (Yusef-Zadeh et 
al]|2006| [Li et al. 2009 ) . 



The first report of periodic behaviour in SgrA* radio 



emission was made by Falcke (1999), who found a 57 day 
periodicity analysing the structure function of the 2.3 GHz 
light curve obtained with the Green Bank Interferometer 
(GBI). At that time SgrA* was also monitored at 8.3 GHz, 
but no periodicity was found at this frequency. At smaller 
wavelengths, analyses of Very Large Array (VLA) observa- 
tions at 3.6, 2.0 and 1.3 cm lasting 20 years, detected a 106 
day cycle in the SgrA* emission (Zhao et al. 2001). In this 
work the period was determined using the Fourier Method 
and the mean profile of the 106 day cycle was obtained af- 
ter 13 cycles. The periodic behaviour was more evident at 
1.3 cm, and there was an increase of variability amplitude 
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with frequency. The 57 and 106 days periods are resonant 
but after the first detection, no other observation was able 
to confirm them, at these or other radio frequencies. |Her-| 
rnstein et al. (2004) also used the VLA to study the long 
term radio variability of SgrA* at 2.0, 1.3 and 0.7 cm on a 
3.3 year campaign. They confirmed the increase of variabil- 
ity amplitude with frequency, but the observations did not 
have sufficient temporal coverage to detect periodicity. 

Several models have been invoked to explain the spec- 
tral energy distribution of SgrA*, as the quasi-Bondi-Hoyle 
accretion f Melia|1992||1994||Coker et al 11999] [Moscibrodzka 
et al. 2006i Cuadra et al.|2006|), disks ([Narayan et al.|1998| 
Yuan et al.[|2003t and jets ( [Falcke Markoff|2000| [Markoff 
et al. 2007|). Two models must be mentioned that try to 



explain the short timescale variability: the Plasmon model 
( van der Laan|1966| |Yusef-Zadeh et al.|2006 r2008) and the 
Hot Spot model (iBroderick k Loeb||2005| 12006), while the 



longer timescale variability can be modelled using the jet 
and/or disk models. As only one model was not sufficient 
to explain all the observational characteristics of the SgrA* 
emission, symbiotic models were proposed, as the jet-disk 
models ([Yuan et al.|2002 ; Kuncicfc Bicknell|2004[|2007|[jor 



2005). Macquart & Bower (2006) even proposed electron 



ley k, Kuncic||20()8 ) and the halo model (Prescher & Melia 



scattering as an extrinsic origin to explain part of the vari- 
ability at radio wavelengths. 

In an attempt to understand the emission of the inner- 
most accretion ffow, recent VLBI observations at 1.3 mm 



inferred 37/ias for the size of the black hole shadow ( |Doele- 
man et al.|2008| . VLBI observations at shorter wavelengths 
do not suffer the inffuence of interstellar scattering and could 
be used to detected the event horizon and the structure of its 
surrounding region ( Falcke et al.||2000 ). However, the emis- 
sion mechanisms and the variability at centimetre and mil- 
limetre wavelengths will only be understood through moni- 
toring programs. 

In this paper we present single dish monitoring of the 
7 mm emission of SgrA* during the last six years. In Sec- 
tion [2] we describe the observations and the method of data 
reduction, in Section |3] we show the light curve and discuss 
variability and periodicity . Finally we present our conclu- 
sion in Section m 



2 OBSERVATIONS 

The 7 mm observations were performed between 2006 and 
2012 with the 13.7 m radome enclosed Itapetinga radio tele- 
scope, localized in Sao Paulo, Brazil. At this wavelength, 
the antenna half power beam width (HPBW) is 2.2 arc min 
and the radome transmission is 0.68. The receiver, a room 
temperature K-band mixer, has a 1-GHz double side band, 
with a ^ 700 K of noise temperature. The calibration was 
made with a noise source of known temperature and a room 
temperature load, which take into account the atmospheric 
and radome attenuation ( Abraham k, Kokubun|199"2 ). Virgo 
A, which has a ffux density of 11 Jy at 7 mm, was used as 
the primary ffux calibrator. 

The method of observation was scans in right ascen- 
sion, with 30 arcmin amplitude and 20 s duration. Only 
24 arcmin were effectively utilized, the rest was lost dur- 
ing the inversion of the scan direction. During each scan 
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Figure 1. Map of the Galactic centre (5^ x 5^) at 6 cm ( |Yusef-| 
[Zadeh et al.|2000| , indicating the the beam size of the Itapetinga 
radio telescope (red circle) and the scan direction (pointed red 
lines). The position of SgrA* is indicate by a black cross. 



81 equally spaced data values, with integration time of 1/6 
s, were obtained. Each observation consisted of the average 
of 30 scans, preceded by a calibration; during an observing 
day, between 7 and 14 observations were obtained and av- 
eraged, always when the source was in elevation above 30° 
and below 75°. The observations of SgrA* were alternated 
with similar observation of the HII region SgrB2, which was 
used as a quasi-simultaneous secondary calibrator. 

The main problem in the detection of radio variability 
from SgrA* with single dish observation is the contribution 
of the continuum emission from the region surrounding this 
point source. As we show in Fig.[l] one scan with 24 arcmin 
amplitude and 2.2 arc min HPBW, detects not only SgrA* 
but part of the continuum complex SgrA. Since only SgrA* 
emission is variable, we used SgrA as an instantaneous cal- 
ibrator, and as a result increased the variability detection 
capacity of our observations. To check the accuracy of this 
calibration and obtain a realistic value of the errors involved 
in the ffux determination, we used a second calibrator: the 
strong HII region SgrB2 Main. 

In order to separate the contribution of SgrA* from the 
surrounding SgrA emission, we ffrst averaged all the obser- 
vations from the year 2007, and adjusted three gaussians to 
the data: a point source with HPBW equal to the beam size 
(2.2 arc min), with two free parameters (peak ffux density 
and central position) representing the average ffux density 
of SgrA* during this period, and two other gaussians with 
completely free parameters. The result of the best fftting is 
presented in Table 1; in Fig. [2] we show the averaged data 
(dots), the adjusted gaussians and the residuals (observa- 
tions minus model). The blue and the navy coloured regions 
in Fig.[l]are represented in Fig.|2]by the two strongest gaus- 
sians; the adjusted HPBWs match very well their angular 
sizes. The emission of SgrA* is represented by the smaller 
gaussian. 
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Figure 2. Averaged 7 mm observations (points) during 2007. 
The pointed and small broken lines represent the gaussian model 
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for the SgrA complex while the large broken line represent the 
Gaussian contribution of SgrA*. The full line is the sum of the 
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Table 1. Parameters the Gaussian models of SgrA*, 
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During the data reduction process, we also adjusted 
three gaussians to the observations, with the HPBW of each 
of them given in Table 1, maintaining constant their rela- 
tive positions as well as the relative intensities of the two 
gaussians that represent the contribution of the SgrA com- 
plex. The free parameters were the position and intensity of 
the strongest gaussian, and the intensity of SgrA*. The flux 
density of SgrA*, calibrated with the SgrA complex, was 
obtained by dividing its observed flux density by the ratio 
of the flux densities of the strongest component of the SgrA 
complex and of its 2007 average value. 

The observations of SgrB2 Main, the quasi- 
instantaneous calibrator, were made with the same 
method as those of SgrA*, i.e. scans in right ascension with 
30 arcmin amplitude and 20 s duration. The observations 
were made alternately with those of SgrA* and with the 
same integration time. As for SgrA, two gaussians were 
fitted to the averaged observations obtained in 2007, their 
parameters are presented in the last two rows of Table [2] 

The flux density of SgrA* was also calibrated using the 
strongest component of SgrB2 Main and he difference of of 
the results of the two calibrations was taken as the measure- 
ment error. 



3 RESULTS AND DISCUSSION 

Although interferometric observations provide flux density 
measurements with better error bars due to the higher angu- 
lar resolution, single dish light curves can have better tem- 
poral resolution for variability studies, as a consequence of 
more available observation time. This is particularly impor- 
tant for SgrA*, because the previously detected periodicities 
are of monthly timescale and so, they need long time moni- 
toring to be conflrmed. Also, the large day-to-day variability 
at 22, 43 and 86 GHz, recently reported by Lu et al. (2011 ), 
needs a sequence of observations in consecutive days to be 
detected. This is the flrst time that SgrA* variability is stud- 
ied using single dish observations and again, this was only 
possible because we used the emission of the SgrA complex 
as a simultaneous calibrator, as discussed in section [2] In 
Table 2 we show the values of the flux density for the flrst 
year of our monitoring. The values for the complete light 
curve are presented in the online version. 



3.1 Variability 

The results of six years monitoring of SgrA* at 7 mm are 
presented in Fig. |3] We detected an increase in the emission 
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Figure 3. 7 mm light curve of SgrA* 
SgrB2 main 



The error bars represent the differences obtained using the two cahbrators, SgrA complex and 



after the first semester of 2009; between 2007-2009 the av- 
erage emission was 2.1 ± 0.2 Jy, while it was 2.8 di 0.2 Jy in 
2010 and 3.1 di 0.2 Jy in 2011. At the same time there was 
an increase in the amplitude of the variability, which had an 
average value of 1.1 Jy in 2007-2009 (low phase) and 2.1 Jy 
afterwards (high phase). The maximum flux density in our 
light curve was 4 Jy, value never reached at 7 mm, although 



Yusef-Zadeh et al. (2011) detected 3 Jy with the VLA before 



the beginning of our observations. 

Coincidence in flux density during quasi-simultaneous 
observations with interferometric techniques, at the same 
frequency, shows the reliability of our results. In fact, |YuseF] 
Zadeh et al. ( 2006 ) , using the VLA for rapid variability stud- 



ies, reported a 43 GHz flux density of 1.70 =b 0.10 Jy on 
2006 July 17, whi le we found 1.90 i fc 0.01 Jy on July 18. At 
other frequencies. Fish et al. ( 2011^ found a flux density of 
2.07±0.15 Jy at 1.3 cm on 2009 April 5, while for this day we 
found 2.21 zbO.13 Jy at 7 mm. Finally, we have 4 days data 
coincident with ATCA 3 mm observations ( Li et al^|2009 ), 
which show the same variability pattern, as can be seen in 
Fig. [4] This is the first time that correlation between vari- 
ability at these two radio frequencies is detected in SgrA*. 
The ATCA observations also revealed the existence of IDVs 
on 2006 August 13. When we looked at our data in a 10 min 
timescale, we found during the same day a dispersion higher 
than average; unfortunately, the large error bars of the ob- 
servations for short integration times makes it difficult to 
detected these events without simultaneous interferometric 
observations. 

In our single dish light curve we have 183 observations 
in consecutive days. We found day-to-day flux density varia- 
tions larger than 0.5 Jy in 23 of them, only four times during 
the low state. In 6 opportunities, the day-to-day variations 
were even larger than 0.75 Jy. Similar variability was de- 
tected in 2007 at 43 GHz using VLBI observations, when 
SgrA* emission showed changes in flux density of 0.5 Jy 
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Figure 4. ATCA and Itapetinga simultaneous observations of 
SgrA*. 



that coincided with two NIR flares ( Kunneriath et al.|2010 



Lu et al. 2011). Unfortunately, there are no infrared data 



for the epochs at which we detected strong daily variability. 
On the other hand, there are 45 consecutive days in which 
the flux density variation was smaller than our error bars. 
These events were present both before and after the rise in 
the light curve, being 16 detected before 2009 December. 

In our monitoring program, we have 5 epochs with more 
than six consecutive days of observation, they are shown in 
Fig. [5] Day-to-day variability is present in some of them; 
the largest occurred in 2010 January, when the flux density 
changed from a minimum of 1.9 =b 0.2 Jy in January 29 to a 
maximum of 3.7 ±0.3 Jy after 4 days. This happened at the 
beginning of the high phase in the light curve. Four months 
later, in 2010 June, the flux density slowly decreased from 
3.2 ± 0.2 Jy to 2.6 ± 0.1 Jy during days 6 and 12. 

The increase in the 7 mm flux density during the last 
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Figure 5. Light curve of SgrA* for five epochs with observations 
in more than six consecutive days. 



years was also observed at other wavelengths. The RXTE 
(1.5-12) keV ASM light curve (Rossi X-Ray Time Explorer, 
All Sky MonitoiQ shows an increase in the SgrA* flux den- 
sity by a factor of 2.5 relative to its average value before 
2010. SgrA* was never so intense at X-ray since the be- 
ginning of the ASM program. In the infrared, the K band 
intensity also showed an increase that started in 2008, af- 
ter an intense flare detected in August 5 by the Very Large 
Telescope ( Dodds-Eden et al.|2011 ). 

We were able to determine the spectral index a {Su oc 
jy^) between the 7 mm and H-band flux densities in two 
simultaneous observations. In 2006 June 24 and 27, the H- 



band flux densities were 4 and 7 m Jy, respectively, ( Bremer 



et al. 2011 ) while the 7 mm flux densities were 1.96=b0.01 and 
2.45 =b 0.05 Jy respectively; between those days there was a 
monotonous rise in flux at 7mm. The IR-radio spectral index 
was -0.7 for both coincident days, which is compatible with 
the IR spectral index of —0.6 ±0.2 found by |Hornstein et ah] 
( 2007 ) . This coincidence favours the idea that the emission 



could be produced by the the same non-thermal electron 
distribution (Yuan et al. 2003). We have also 7 mm and 



K-band simultaneous observations on 2006 June 21, which 
gave a spectral index of -1.1. 



3.2 Periodicity 

As it was already mentioned, two periods were previously 
reported in the SgrA* radio emission: 106 days, determined 



from 20 years of VLA observations at 22 GHz ( Zhao et al. 
2001 ), and 57 days found after 600 days of observations with 
the GBI at 2.3 GHz ( |Falcke|[l999t . Our data set is not so 
long as the VLA database, but has better temporal resolu- 
tion, while it is longer than the GBI monitoring but with 
worse temporal resolution. In other words, our light curve 
is irregularly sampled but without long gaps, which allows 
a more reliable statistical analysis. 

To check the existence of possible periodicities we 
used the Structure Function Analysis and the Stellingwerf 
method. The intensity of the structure function (hereafter 
SF) is deflned as D{r) = {[I{t + r) - /(t)]^). This is a sim- 
ple method to verify how much the intensity / varies after 
a time delay r. If there is no characteristic timescale in the 
light curve, D will have a random value for each r, and so 
there will be no maxima or minima in the structure function. 
Instead, if there is a typical timescale in the light curve, 
D will be large when r — and appear as a maximum in 
the SF. This method is also useful to flnd periodicity, be- 
cause if there is a period P, D will have a low value when r 
is equal to an integer number of periods, and this will result 
in minima in the SF. 

The Stellingwerf method (Stellingwer f||1978t is used to 
flnd periodicity in an irregularly sampled series. The test 
divides the sample in m groups following the phase criterion: 



U 
P 



(1) 



where 02 is the phase, U is the observation time and P is 
the guessed period. The brackets correspond to the integer 
part of the fraction. Each observation day of the sample has 



http : // xte.mit.edu/ ASM _lc.html 



6 P.P.B.Beaklini and Z. Abraham 



a value for the phase and it is grouped with other days with 
similar phases. The variable is defined as the ratio be- 
tween the sum of the mean square deviation of the fluxes in 
each group, and the mean square deviation of the complete 
sample. If P did not correspond to the real period, will 
be close to 1. If P is indeed the period in the signal, there 
will be in the same group days with similar flux densities, 
and so, the sum of the group mean square deviation will be 
smaller than that of the complete sample. The number of 
groups is a free parameter, but higher values of m provide 
deeper minima in the 6 vs. P curve. However, m can not be 
too large, because we need a sufficient number of elements 
in each group to obtain a reliable statistics. We used m = 30 
for all tests in this work. 

First, we applied both tests to the original light curve; 
the results are shown in Fig.|6](a) and (b). Then, we applied 
the tests to the data after removing a third-order polynomial 
to discount the slow time variation; the results shown in Fig. 
6|(c) and (d ). This last procedure is similar to that used by 
Zhao et al.l (12001). 



In both cases we found several maxima and minima. 
To verify their significance and check if they are not con- 
sequence of sampling problems, we constructed an artificial 
light curve using random values for the flux density but with 
the same variability amplitude, for the days in which SgrA* 
was observed. The results of the statistical tests applied to 
this light curve are presented in Fig. [7| The range of the os- 
cillations in the SF function lies between —0.23 and +0.43 
of its mean value, which corresponds to 2.5 and 4.0 a, re- 
spectively. We used these values to limit the significance of 
the maxima and minima in Fig. [6](c), shown as red dashed 
lines. The only minimum deeper than this limit corresponds 
to 156 it 10 days, but marginal minima are present at the 
resonant 320 ± 10 days period, and at 220 it 10 days, which 
could be resonant with the 106 day cycle reported by |Zhao] 



et al. (2001). 



The significance of the Stellingwerf method results can 
be computed through the F-Test, developed by |Kidger et al.| 
(T992I: 



e 



e 



(2) 



where / > 0.5 indicates strong periodicity and / < 0.25, the 
absence or weak periodicity. In fact, the / test applied to 
the random light curve results shows always / < 0.15. This 
value is represented as the red dashed line in Fig.[6](d). We 
can see wide 156 and 320 days minima which seem to be 
significant (/ = 0.27), as well as several other very narrow 
minima, which we interpret as fluctuations. As discussed in 



Kidger et al. (1992), the Stellingwerf method can not reveal 
periods in a sampling smaller than six times the proposed 
period. In our light curve, this means periods longer than 
350 days. 

We conclude from the Structure Function and the Stel- 
ingwerf method that are two possible periods: 156 it 10 and 
220 ± 10 days that do not coincide with any of those previ- 
ously reported: 57 days ( Falckep999 ) at 2.3 GHz and and 
106 days (Zh ao et al.|2001| ) at 22 GHz. However, we specu- 
lated if the existence of two beating periods could be the rea- 
son for the different results obtained by different authors. To 
check this hypothesis, we constructed artificial light curves 
as combination of two sine square functions, leaving their 



relative amplitudes as a free parameter and adding a ran- 
dom signal with amplitude similar to that of the periodic 
functions. We tried all possible combinations of 57, 106, 156, 
212, 220 periods. We used 212 days because it is exactly a 
resonance of the 106 day cycle. Two combinations repro- 
duced our statistical results: 156 and 212 days, and the 57 
and 156 days; the amplitudes of the 156 days periodic sine 
square functions were 1.2 and 1.5 times larger than that of 
the 57 and 212 days periodic signals, respectively. The sta- 
tistical tests are compared in Fig.js] The 156 days minimum 
was not found by any combination of other periods, and the 
result of using 220 or 212 days period are very similar. 

The minimum that would correspond to the 57 days 
period did not appear in any of the combinations, even if this 
period was used in the input data. We believe that this could 
be a consequence of not having enough temporal resolution 
added to the existence of large amplitude short timescale 
fluctuations. Therefore, we repeated our calculations using 
our simulated light curve for the 57 and 156 days periods 
with daily data points, and in fact, the presence of the 57 
days minimum can be seen in Fig. [9] It is important to 
notice that the SF of the simulated light curve is in very 
good agreement with the SF found by [Falcke (1999) for the 
2.3 GHz data and that the larger amplitude of the second 
peak seen in both SFs is a consequence of the existence 
of these two periodicities. This was not the case when the 
combination of 156 and 212 days periodic signals were used, 
even with daily resolution. 

In other words, only the combination of 57 and 156 days 
reproduce both our statistical tests and Falcke's. To verify 
the range of possible periods around 57 days that can still 
reproduce Fig. [S] we repeated the tests combining the 156 
day period with different periods between 51 and 61 days. If 
the second period is smaller than 55, there is no minimum 
in the SF around 220 days. Otherwise, if the second period 
is larger than 59 days, the minimum in 156 days is not deep 
enough. This indicates that a periodic signal of of 57it2 days 
mixed with other of 156 days can reproduce the 156±10 and 
220 dt 10 minima revealed by our statistical tests. 

Considering the existence of the 57 it 2 period, we can 



interpret the 106 day cycle reported by Zhao et al. (2001) 



as a resonance of that period. To check this hypothesis we 
modulated our light curve with a 106 day period, and com- 



pared the result with that of Zhao et al. (2001). We found 



a similar pattern for both cycles, but without the two large 
peaks on days 75 and 100 present in the 13 mm cycle. The 
peaks can be due to the low number of observation days in 
the VLA light curve. 

Finally, we must mention that no periodicity was re- 
ported by Macquart Bower (2006) for the VLA data ob- 



Herrnstein et al. (2004). To understand that we 



tained by 

calculated the SF for these data and we did not find any 
periodicity either. Again, the reason is that the VLA light 
curve has not enough time coverage to detect monthly peri- 
odicities. 



3.3 The infalling ionized gas cloud 

As presented in subsection 3.1, variability of SgrA* was de- 
tected at 7 mm in timescales of days and years. The lat- 
ter consisted of a continuous increase in flux density that 
started in 2008, as can be seen in FigjlT] where the radio 
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Figure 6. Statistical results for SgrA* 7 mm light curve: (a) and (b) SF and Stellingwerf results for the original light curve; (c) SF 
and (d) Stellingwerf and F-test results (in red) after a third-order polynomial was removed from the light curve. The red dashed lines 
represent the significance limit. 
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Figure 7. Statistical results for a random signal in the same days as those of our light curve. For the SF, the red dashed lines represent 
the extremes of the oscillations; for the Stellingwerf method the F-test function is included (in red), the dashed red line is the maximum 
value of this function 



data are averaged in 60 days bins. Coincidently, [Gillessen et| 
al. (2012) detected a dense gas cloud that is moving in the 



direction of SgrA* in a highly eccentric orbit and will reach 
its minimum approach in 2013.5. The almost spherical cloud 
of 15 mas (2 x 10^^ cm) observed in Br7 in 2008, is being 
tidally disrupted and in 2011 presented a tail elongated in 
the orbital direction of about 200 mas. 

To investigate if the increase in the 7 mm flux density 
between 2008 and 2011 could be the result of the appearance 
of the disrupted tail, we calculated the flux density of an 
ionized optically thick cloud of temperature T that subtends 
a solid angle — ttO^: 



Sx 



2kT 



(3) 



where Sx is the flux density at wavelength A, and k is the 
Boltzmann constant. 

Assuming T = 10^ K and 6> = 150 mas, the size of the 
observed tail, we obtain Srmm — 1 Jy, in agreement with 
our observations. We could assume then that the increase in 
the observed 7 mm flux density between 2008 and 2011 cor- 
responded to an increase in the solid angle Q. This increase 
seemed to stop around 2011, which could be interpreted as 
the cloud becoming optically thin at that epoch. Assuming 
the optical depth rrmm = 1, we can calculate the emission 
measure EM — nj^L of the cloud from: 



TA = 2 X 10 rinLj^Qff 



(4) 



For T = 10 K, the Gaunt factor is ^// ^ 7, and EM 
1.4 X 10^^ cm-^ 



If we assume a density of 10^ cm~^ for the cloud, as 
obtained from numerical simulations of the evolution of a 
disrupted cloud ( Burkert et al.||2012 ), we obtained a length 
L = 1.4 X 10^"^ along the line of sight, showing that the tail 
in the orbital plane could be formed by a thin layer of gas. 

If the 7 mm radio emission is due to free-free emission, 
we should expect a corresponding flux density in the IR, 
given by 

Sx = 5AxlO-'^T-°'gffi\,T)QEMeM^) Jy- (5) 

For the IR, gff{X,T) ^ 1; resulting for the Ks band 
a flux density of Sks ^ 0.057 Jy. Although this value is 
much higher than any flux density detected in the SgrA* 
region, it must be noticed that it corresponds to an extended 
emission, therefore the signal should be diluted in an area 
of 150 X 150 mas^, which represents the size of the tail. 
Considering a pixel size of 13 or 27 mas per pixel, typical 
of the VLT data (Dodds-Eden et al. 



2011 



Gillessen et al. 



2012| ) we obtain 0.43 and 1.8 mJy per pixel, which could 
have been overlooked in the data. 

However, the increase in the amplitude of the fluctua- 
tions in timescales of days simultaneously with the increase 
in the total flux density that started around 2009, favours 
the hypothesis that is intrinsic to SgrA*. In fact, this relation 
between rms fluctuations and flux density was already de- 
tected in active galaxies and X-ray binaries a t X-rays (^Utt-| 
ley McHardy|2001| )and optical wavelengths ( |Gandhi|2009f 
when they are in a luminous phase, which lead |Heil et al.| 
(2012) to claim that this can be an overall propriety of the 
accretion flow in compact sources. In the case of SgrA*, the 
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Figure 8. Comparison between the statistic results for the real (upper) and simulated (lower) light curves, for two different period 
combinations: 156 and 212 days (left) and 57-156 days (right). Top: Structure Function intensity, bottom: Stellingwerf method. 



increase in the accretion rate could be related to the influ- 
ence of the infalling cloud or to the winds of faint stars dur- 
ing their periastron passage ( Loeb|[2004 ), as those reported 



by Dodds-Eden et al. (2011) 



4 CONCLUSION 

In this paper we presented the results of single dish mon- 
itoring of SgrA* at 7 mm during a six year period. Great 
care was taken to separate the contribution of SgrA* from 
the surrounding emission, the latter used, together with the 
HII region SgrB2, as an instantaneous calibrator. The reli- 
ability of our results was verified by comparing them with 
interferometric observations at coincident epochs. 

During the observing campaign, several epochs with 
consecutive day observations allowed the study of variabil- 
ity in daily as well as longer timescales. The resulting light 
curve shows a slow variation in a timescale of years, a ran- 
dom component with shorter timescale and a day-to-day 
variability with large amplitude variations between two con- 
secutive days. 

We summarize the variability result as follow: 

(i) SgrA* emission at 7mm has been growing up through 
the years since 2008, this behaviour coincides with an 
increase in X-ray emission, and there is also some evidence 
of high activity in the infrared stating at the beginning of 
2009. 

(ii) SgrA* light curve presented episodes of high 
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Figure 9. Top: SF results for the simulated two period light 
curve (57 and 156 days) with daily time resolution. Bottom: the 
SF results obtained by 
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Figure 10. 106 day modulation of the 7 and the 13 mm hght 
curves 
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Figure 11. 7mm hght curve with 60 days binning and the third 
order polynomial removed from the light curve (continue line) 



amplitude day-to-day variability that were more common 
in the high phase; 

(iii) The variability at 7 mm showed the same pattern 
than that at 3 mm during four coincident days in 2006; 

(iv) Simultaneous 7 mm and IR data revels that the 
radio-IR spectral index is similar to the IR spectral index, 
suggesting that the emission at that epoch was probably 
produced by the same electron population; 

We investigate the possibility that the increase at 7 mm 
flux density that started in 2008 could be related to the ap- 
pearance of a Br7 emitting gaseous tail of tidally disrupted 
material, arising from the compact gas cloud that is falling 



towards SgrA* in a highly eccentric orbit |Gillessen et al 
|2012 ). Even though it is possible to explain the 7 mm flux 
density increase as free- free emission from such a cloud, the 
predicted IR emission was not observed. However, the in- 
crease of the fluctuations together with the flux density at 
7 mm , favours the intrinsic origin of this extra emission 



and is compatible with the rms/flux relation found in other 
compact sources ( [Uttley k, McHardy|2001 ). 

In search for periodicity, the light curve was statistically 
analysed using the structure function and the Stellingwerf 
method. The result can be summarized as follows: 

(i) The statistical tests realized in this work revealed the 
existence of two minima: 156 di 10 and 220 it 10 days, but 
did not show any minima related to the 57 or 106 day 



periods reported by Falcke ( 1999 ) and Zhao et al. 



(|200l').To test the reliability of our results we constructed 
an artiflcial light curve with random signals in the same 
observation days; no periodicity was found, showing that 
the existence of periodicity was not due to sampling effects. 

(ii) We investigate if the divergence between the detected 
periods can be a consequence of the existence of two 
periodic signals. For that we simulated a light curve, as a 
combination of two quadratic sinusoidal functions, with 
data at the same observation days as those of our light 
curve, and verifled which periods and relative amplitudes 
better reproduced our statistical results. We concluded 
that two combinations of periods can explain our 
statistical results: 156 plus 212 days, and 57 plus 156 days. 
The amplitudes of the 156 days periodic functions were 1.2 
and 1.5 times larger than that of the 57 and 212 days 
periodic signals, respectively. 

(iii) To test the influence of sampling in the detected 
periodicity, we constructed a simulated light curve for the 
57 and 156 days periods with daily data points, using the 
same relative amplitudes as in the previous best test. The 
result of the structure function was in very good agreement 
with that found by [Falcke | ( [l999| for 2.3 GHz, which 
revealed the 57 day periodicity. 

(iv) We also modulated our light curve with a 106 day 
period, to compare it with [Zhao et al. (2001) results at 13 
mm. The modulated light curves shows a similar pattern 
but without the two large peaks on days 75 and 100 of the 
cycle. This similarity can be explained as a consequence of 
a resonance of the 57 day period. 

We conclude that two periodic signals added to a ran- 
dom component can reproduce statistical results similar to 
what it was obtained for the SgrA* emission. The combina- 
tion of 156 ± 10 and 57 =b 2 day periodicity is what better 
explains our results. If real, these two periods may be ex- 
plained by precession and/or nutation of a accretion disk. 
The precessing disk in SgrA* could be consequence of the 
the Bardeen- Peterson effect ([Bardeen k, Petterson|1975 ), as 
discussed by Liu & Melia (2002). In fact, this process can 
lead to a hundred day period in the SgrA* radio emission 



( .Caproni et al., 2004, , Rockefeller et al.|2005 ) 
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